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Based on a localized crystal electric field model for the U 4+ in the (5/) 2 -configuration, we 
analyze the resonant x-ray scattering spectra around U Miv and My edges in URu2Si2, taking 
full Coulomb and spin-orbit interactions into account. We consider two level schemes, a singlet 
model of Santini and Amoretti and a doublet model of Ohkawa and Shimizu, and assume 
the antiferroquadrupolar order and the antiferromagnetic order as candidates for the ambient 
pressure phase and the high pressure phase. It is found that the spectral shapes as a function 
of photon energy are independent of the assumed level scheme, but are quite different between 
the antiferroquadrupole and antiferromagnetic phases, This may be useful to determine the 
character of the ordered phase. 
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1. Introduction 

The ternary intermetallic compound URu2Si2 has 
been attracting much attentions since the discovery of 
the coexistence of antiferromagnetic (AFM) ordering be- 
low To = 17.5 K and superconductivity below Tq = 1.2 
K. 1 The ordering pattern of the AFM order is type-I 
structure with the propagation vector Q = (001) along 
the c-axis (Fig. 1). The phase transition at To is charac- 
terized by bulk anomalies such as specific heat, 2 linear 
susceptibility, 2 non-linear susceptibility, 3,4 thermal ex- 
pansion 5 and electrical resistivity, 6 strongly indicating 
that the transition is of second order. 

The observed magnetic moment is unusually tiny 
(Mord ~ (0.03 ± 0.01)/xb per U ion at saturation). 7-9 It 
is too small to account for the observed bulk anomalies. 
The magnetic excitation spectra shows energy gap be- 
low To according to neutron scattering experiments. 10 
Under the magnetic field, the magnetic moment and the 
magnetic excitation gap exhibit the different field depen- 
dencies. 11 Based upon these observations, it is widely 
believed that, the AFM order below To is merely a sec- 
ondary order and there exists some unknown primary 
order, so-called "hidden order (HO)". 12 As for the true 
nature of the HO phase, various competing theoretical 
possibilities have been demonstrated. For instance, the 
antiferroquadrupole (AFQ) order, 13-17 the multispin cor- 
relation, 18 the unconventional spin density wave, 19 inho- 
mogeneous orbital current, 20-22 Jahn- Teller distortions 23 
and so on were proposed. However, it is safe to say that 
the microscopic mechanism governing the whole phenom- 
ena is still in a matter of controversy. 

Recently, neutron scattering and 29 Si NMR experi- 
ments were performed under pressure. 24-29 The mea- 
surements are interpreted that not the magnitude of the 
AFM moment but the volume fraction of the AFM phase 
is tiny below To. The volume fraction of the AFM phase 
is estimated about 1 % at ambient pressure while the 
magnitude of the moment is 0.25 /zb per U ion. Another 
discovery of importance is the presence of the phase tran- 
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Fig. 1. Crystal structure of body-centered tetragonal URu2Si2- 
The arrows show the magnetic structure corresponding to a wave 
vector Q = (001). 

sition around pc = 1.5 GPa, accompanied by the discon- 
tinuous changes of the quantities such as To, the lat- 
tice constant and the value of the magnetic moment. In 
high pressure phase, the AFM order extends uniformly 
throughout the entire sample. It exhibits a three dimen- 
sional AFM Ising behaviors with the saturated staggered 
moment ~ 0.4/ia/U and with the same magnetic struc- 
ture as in the ambient pressure (HO) phase. 24 

One of the promising efforts to elucidate the nature 
of the order parameter in the d and / electron systems 
is the resonant x-ray scattering (RXS) measurement. By 
utilizing the facts that the method is element specific 
and level specific, and that the RXS intensity generally 
shows special incident photon polarization and azimuthal 
angle dependencies, we can deduce valuable information 
to understand the physical properties of subjects. In 4/ 
systems, for example, the RXS peaks at the Ce Tin ab- 
sorption edge in CeBg are interpreted as the direct obser- 
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vation of the AFQ ordering 30-33 while the peaks at the 
Dy Liu edge in DyB 2 C2 are explained as the consequence 
of lattice distortion. 34-37 

For actinide compounds, the Mpv.v edge resonance is 
available for supcrlattice Bragg spots. Since the dipolc 
process gives rise to a transition from the 3d states to 
the 5/ states which drive ordering, the RXS signal is 
more direct than those at the K-edge in the transition- 
metal compounds. In transition metal oxides, the RXS 
signals at the K-edge involve the transition from the Is 
state to the Ap states. Since the 4p states are considerably 
extending in space, they are very sensitive to electronic 
structures at neighboring sites. Thereby the RXS signal 
on the orbital ordering superlattice spots are mainly con- 
trolled by lattice distortion. 38-40 

Isaacs et al. reported the resonant enhancements cor- 
responding to the U M\y edge upon URu2Si2 at ambient 
pressure below T . 41 The measured RXS peak intensity 
at the prohibited Bragg spot (003) develops linearly over 
an unusually wide range of temperature from T down 
to about 3 K, then ceases to grow when superconductiv- 
ity sets in. The width of the spectra is evaluated about 
5 eV. Later, Lidstrom et al. have observed the peak in- 
tensities at U Miv and My edges from Ui- x Np a; Ru2Si2 
alloys with x = 0.1,0.5 and 1.0 also using the RXMS 
technique. 42 

In this paper, we analyze the RXS spectra in the low 
temperature phases of URu2Si2 on the basis of a localized 
model under the crystal electric field (CEF). Although 
the CEF-like excitations are not clear in the high tem- 
perature phase, a localized picture can explain well the 
behavior of various quantities such as the specific heat 
and the magnetic susceptibility. 13 ' 24,25 Since the CEF 
parameters are not known, the level schemes are not def- 
initely determined. There are two prevailing schemes, the 
"singlet" scheme of Santini and Amoretti 13 and the "dou- 
blet" scheme of Ohkawa and Shimizu. 17 On the basis of 
the two level schemes, we consider both the AFQ and 
AFM phases, and calculate the RXS spectra as a func- 
tion of photon energy and of azimuthal angle. We find 
that the calculated spectral shapes are independent of 
the assumed scheme, although the intensities are quite 
different. 43 Thereby the RXS spectral shapes depend 
only on the assumed order parameter. The RXS spectral 
shapes as a function of energy are found quite different 
between the AFQ phase and the AFM phase. In princi- 
ple, this may make it possible to distinguish whether the 
experimental signals come from the HO phase or from 
the AFM domains. Unfortunately available experimental 
spectra are not clear enough to distinguish those differ- 
ences. 41,42 Shishidou et al. 44 calculated the RXS spec- 
tra in the AFQ phase within the doublet scheme. The 
present calculation agrees with their result. 

The present paper is organized as follows. In §2, we 
briefly summarize the theoretical framework of RXS. In 
§3, we discuss two level schemes under the CEF on the 
basis of the localized 5/ model, and the AFQ and AFM 
states as candidates of the ground state under these 
schemes. The intermediate states are briefly discussed. In 
§4, numerical results of the RXS spectra are presented in 
comparison with the experiments. A summary and dis- 




Fig. 2. Geometry of the RXS experiment. Photon with polariza- 
tion c or 7r is scattered into the state of polarization <r' or n' at 
Bragg angle 9. 



cussion are given in §5. In Appendix, general expressions 
of the matrix elements of the scattering amplitude are 
given. It is proved that the RXS spectral shapes are in- 
dependent of the assumed level schemes. 

2. Theoretical Framework of RXS 

We consider the experimental situation of RXS shown 
in Fig. 2. The incident photon with frequency w and wave 
vector ki is scattered into the state with frequency lo 
and wave vector kj. The sample is rotated around the 
scattering vector G(= k/— lq) by azimuthal angle ip. The 
resonant enhancement is found at the U M IV and My 
absorption edges, where the 3d core electron is virtually 
excited to the 5/ states then falls to the core level in the 
dipolar (E\) process. The tuning energy for M\y edge 
is around 3.728 keV and that for My is around 3.550 
keV. The scattering tensor can be well approximated by 
a sum of the contributions from each site of the created 
core hole. Therefore the RXS intensity in the E\ process 
is given by 



(2.1) 



with 



(y-oMj)|A)(A|aw(j)|y> ) 
hu)-{E A - Eq) + iT 



(2.2) 



where N is the number of U ions. The \ip ) represents 
the ground state with energy E , and |A) represents the 
intermediate state with energy E A . The dipole operators 
x a (j)'s are defined as x\{j) = x, x 2 (j) = y, and Xs(j) — 
z in the coordinate frame fixed to the crystal axes with 
the origin located at the center of site j. The P M and P' M 
denote the geometrical factors for the incident photon 
and the scattered photon, whose explicit expressions are 
given in ref. 37. The T describes the life-time broadening 
width of the core hole. 
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Tabic I. Slater integrals and spin-orbit interaction parameters 
within the HFA (in units of eV). 45 

the (3d) 10 (5/) 2 configuration 



Table II. Eigenstates of the CEF Hamiltonian eq. (3.3) with \M) 
denotes the state of J z = M. 
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3. Electronic Structure 

3.1 The ground state 

We adopt a localized description that each U ion is in 
the (5/) 2 -configuration. The Hamiltonian of U ions con- 
sists of the intra-atomic Coulomb interaction between 5/ 
electrons and the spin-orbit interaction (SOI) of 5/ elec- 
trons, which are represented by 91 states in the (5/) 2 - 
configuration. We evaluate the Slater integrals for the 
Coulomb interaction in the 5/ states and the SOI param- 
eters within the Hartree-Fock approximation (HFA). 45 
The values are listed in Table I. It is known that the 
values for the isotropic part of the Slater integral F° are 
considerably reduced due to large screening effects. We 
reduce the value by multiplying a factor 0.25. On the 
other hand, the values for the anisotropic part of the 
Slater integrals are not screened so much. We slightly 
reduce them by multiplying a factor 0.8. We diagonalizc 
the Hamiltonian, and obtain nine degenerate states be- 
longing to the J — 4 subspace as the lowest-energy state. 
Note that these states of J = 4 arc slightly deviated from 
those of the LS coupling scheme with L = 5 and S = 1 
due to the SOI. 

It is known that the crystal structure of URu2Si2 is 
the body-centered tetragonal (bet) ThCr 2 Si2 structure 
(I4/mmm) with lattice constants a = 4.124A and c = 
9.582A at T = 4.2 K (Fig. I). 46 Under the bet symmetry, 
the degeneracy in the J = 4 subspace is lifted by the CEF 
Hamiltonian, which may be expressed as 



Hcry — 



(3.3) 



where O^'s denote the Stevens operator equivalence. 
The CEF parameters B^'s take non-zero values for 
(M)=(2,0), (4,0), (4,4), (6,0) and (6,4). The eigen- 
states of the Hamiltonian split into five singlets and two 
doublets as listed in Table II, where the e, 7, a depend 
on the CEF parameters. At present, it is difficult to de- 
termine the CEF parameters in a convincing way. In the 
following, we discuss two prevailing CEF level schemes. 

3.1.1 Singlet scheme 

Analyzing the temperature dependences of magnetic 
susceptibility and specific heat, Santini and Amoretti 13 
have proposed a CEF level scheme that the lowest three 
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states become r t3 (or r t4 ), T$ and r t2 in order, with en- 
ergy separation 3.8 meV, 9.6 meV from the lowest state, 
respectively, and that the other states have energies more 
than ~ 40 meV higher than the above. This scheme has 
been considered for bringing about the AFQ ordering at 
low temperatures. It successfully reproduced the experi- 
mental magnetic susceptibility and specific heat in a wide 
temperature region. 13 

Within the states T t 3, and T t 2, the dipole operator 
J z and the quadrupole operators O x 2_ y 2 = ( ' - ■ 
O xy = (J x Jy + JyJx) are represented as 




- 7 2N I 
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(3.4) 



(3.5) 



(3.6) 



Another choice of the states T t 4, T t2 makes the ma- 
trix of O x 2_ y 2 and that of O xy interchanged. The e and 7 
are estimated as —0.6775 and 0.2864, according to San- 
tini and Amoretti's analysis. 13 For describing the AFQ 
phase, we add a term XQ(Q) on the basis of the mean- 
field approximation to the multipolar interaction. Here 
(Q) is the staggered moment of one of the quadrupole 
O x 2_ y 2 or O xy and A = 0.185 meV is a mean-field coeffi- 
cient. 13 ' 47 The self-consistent solution gives us the AFQ 
ordering at low temperatures; = 12.1 at zero tem- 

perature. The above expressions indicate that the O x 2_ y 2 
(or O xy ) ordering is most realized at low temperatures 
through the off-diagonal element between the lowest and 
the first excited states. 

For describing the AFM phase, we now introduce a 
molecular field term X' J Z (J Z ), where (J z ) is the staggered 
dipole moment. Since J z has no matrix element to the 
lowest energy level (see eq. (3.4)), the molecular field 
energy has to exceed the crystal field energy for giving 
rise to the AFM phase at low temperature. For A' > 0.468 
meV, a self-consistent solution exists as \(J Z )\ > 3.490 
with h/hb = \(L Z + 2S Z )\ = 2.887. For instance, when 
A' = 0.5 meV, \(J Z )\ — 3.540, and the staggered moment 
becomes fi — 2.925^b- Therefore the AFM ordering can 
be brought about only with assuming a large intersite U- 



4 J. Phys. Soc. Jpn. 



Full Paper 



Tatsuya Nagao and Jun-ichi Igarashi 



U coupling. This staggered moment is much larger than 
the experimental value \i ~ 0.4itB in the pressure range 
between p = 1.5 and 2.8 GPa. 24 - 25 

3.1.2 Doublet scheme 

Ohkawa and Shimizu 17 have proposed another model 
that the lowest energy levels belong to the doublet or 

(2) 

T t5 . This doublet scheme naturally gives rise to the AFQ 
and AFM orderings under moderate intersite U-U cou- 
plings. Also it is consistent with the study of the dilute 
U x Thi_ x Ru 2 Si 2 alloys (0 < x < 0.07) by Amitsuka and 
Sakakibara, which indicates that the lowest CEF state 
of 5/ electrons is magnetically degenerate due to the ab- 
sence of saturation in the susceptibility of U contribution 
down to 100 mK. 48 

Choosing as a lowest doublet, we explicitly write 

|t> = cosa| +3) +sina| - 1), (3.7) 



U M v edge 



U M lv edge 



cos a\ — 3) + sin a\ + 1) , 



(3.8) 



with \M) denoting the state of J z — M. Operators J z , 
O x 2_ y 2 and O xy are represented within the subspace as 



J z = (3 cos 2 a — sin 2 a) 



V% sin a (5 sin a + 3v / 7 cos a) 



O xy — V3 sin a(— 5 sin a + 3\/7 cos a) 



1 
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i 



(3.9) 



(3.10) 



(3.11) 



The intersite interaction induces the AFQ ordering or 
the AFM ordering. 

For the AFQ ordering, the ground state is given by an 
alternating array of the states |+) and |— ), where 



i±>=y|[i r>±u>], 

for the O x 2_ y 2 ordering, and 



l±> 



2H T)±i| I)], 



(3.12) 



(3.13) 



for the O xy ordering. 

For the AFM ordering, the ground state is just an 
alternating array of | |) and | J.). The staggered magnetic 
moment is evaluated from 

fi/fiB - cos 2 a\(+3\L z +2S z \+3)\+sm 2 a\(-l\L z +2S z \-l) 

(3.14) 

Putting a ~ 52.48° ora~ 68.95°, we obtain fi ~ 0.4// B , 
close to the experimental value. 24, 25 

3.2 The intermediate state 

The E\ transition causes a 3d-core electron to an unoc- 
cupied 5/ states, leading to the (3d) 9 (5/) 3 -configuration. 
The Slater integrals and the SOI parameters are evalu- 
ated within the HFA, which are also summarized in Table 
I. 45 The Coulomb interactions and the SOI of 3d and 5/ 
electrons are represented within 364 x(2jd + 1) states 
in the (3(i) 9 (5/) 3 -configuration to construct the Hamil- 
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Fig. 3. Absorption coefficient as a function of the photon energy 
near U Mjv and U My edges. The upper panel is the exper- 
imental result (Ref. 48) and the lower panel is the calculated 
result. 



tonian, and the intermediate state |A) is calculated by 
diagonalizing the Hamiltonian. Here jd = § and | mean 
the total moment of the 3(i-core electron from the Mjy 
and My edges, respectively. The CEF energy is so small 
in comparison with the multiplct energies that it can be 
neglected in the intermediate state. Then the energy lev- 
els are degenerate with respect to the magnetic quantum 
number M within the space of the total angular momen- 
tum J. 

4. Calculated Results 

4-1 Absorption coefficient 

We discuss the absorption coefficient A{u) at U Miy 
and My edges. Within the E\ transition, it is expressed 
as 



(^-£ A ) 2 + r 2 ' 

(4.15) 

We calculate A{u>) by assuming the AFQ ordering or the 
AFM ordering at zero temperature. The final state of 
photoabsorption is the same as the intermediate state 
of the RXS explained in the preceding section. The T 
is assumed to be 0.5 eV. Figure 3 shows the calculated 
result in comparison with the experiment. 49 The core- 
level energy is adjusted such that the calculated peak 
position coincides with that of the experimental peak 
position. The obtained spectra are independent of the 
choice of the ordered phase, which can be verified by 
means of the results obtained in Appendix. The broad 
' peak is due to multiplct structure of the final state. 

4.2 RXS Spectra 
4.2.1 AFQ phase 

Assuming the AFQ phase in the initial state and us- 
ing the intermediate state in the preceding section, we 
calculate the RXS spectra at U M\y and My edges from 
cq. (2.2). The scattering amplitude has a form 
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Fig. 4. The RXS spectra as a function of the photon energy at 
G = (005) with xj) = 0. Left and right panels show the results 
around U My edge and Mjv edge, respectively. Top right panel 
shows the experimental result. 41,50 The rests are calculated re- 
sults corresponding to those in the J z , O x 2_ y 2 and Oxy phases 
from bottom to top, respectively, for both sides of the panels. 



Fig. 5. Azimuthal angle dependence of the RXS peak intensities 
at G = (005) near U Miy edge in the AFQ phases. The solid 
and broken lines represent the results obtained in the Oxy and 
O x 2_ y 2 ordcrings within the doublet scheme, respectively. The 
mixing angle is fixed to being a ~ 52.48° and both a — a' and 
a — 7r' channels arc shown. 



/ 1 \ 

oc 1 , for O xy ordering. 

\ J 

(4.16b) 

Combining these expressions with the geometrical factors 
in eq. (2.1), we obtain the intensities for G = (00£) 

f cos 2 (2^), // = a' 

W«| S i n 2(2^) s in 2 0, (j,' = tt' ' (417) 

for the O x 2_ y 2 ordering, and 

Woc( Sin2 2 K. 2 Q ^ = \ , (4.18) 
M \ cos 2 (2V>)snr 6», M =71" 

for the ordering. 

We find that the spectral shapes as a function of energy 
are independent of the assumed level schemes, although 
the intensities are quite different between the schemes. 43 
Figure 4 shows the RXS spectra at G = (005) with tp = 
near U Miy and My edges, as a function of photon en- 
ergy. For the Miy edge, the spectra have a width of 
~ 5 eV due to the multiplet structure in the interme- 
diate state; a large peak is at fuv = 3.726 keV. and a 
shoulder is at ~ 2 eV higher than the peak. The calcu- 
lated width looks somewhat narrower than the experi- 
mental one. 41,50 For the My, a narrower peak is found 
at hw = 3.549 keV. 

As clear from eqs. (4.17) and (4.18), the RXS intensity 
depends considerably on the azimuthal angle. Figure 5 
demonstrates this dependence. This is in contrast with 
the RXS intensity in the AFM state, where it shows no 
ip dependence as shown below. 

4.2.2 AFM phase 

We also calculate the RXS spectra, assuming the 
AFM phase in the initial state. The scattering amplitude 



M(G,uj) has a form 

/ \ 

M(G,w) oc -£(w) . (4.19) 
\ / 

Combining this with the geometrical factors in eq. (2.1), 
we can verify that the RXS intensities show no signal in 
the a— a' channel while those in the a— ir' channel remain 
constant with varying azimuthal angle ijj. This azimuthal 
angle dependence is valid as long as the scattering vector 
is in the form of (00£). 

As the same as the case of AFQ state, we find that 
the spectral shapes as a function of photon energy are 
independent of the assumed level schemes. 43 The RXS 
spectra at G = (005) near U M\y and My edges are 
shown in Fig. 4. The spectral peak is located at position 
~ 2 eV higher than that in the AFQ state, and the shoul- 
der is seen at the lower energy side. This spectral shape 
is clearly different from that in the AFQ state. Unfortu- 
nately, the experimental spectra are not clear enough to 
distinguish the difference. 

5. Concluding Remarks 

We have analyzed the RXS spectra near U M\y and 
My edges on the basis of a localized CEF model. By 
employing two level schemes, the singlet model of San- 
tini and Amoretti and the doublet model of Ohkawa 
and Shimizu, we have considered the AFQ phase as a 
candidate of the HO phase and the AFM phase as the 
one realized under high pressure p > 1.5 GPa. We have 
calculated the RXS spectra by taking into account the 
full Coulomb and spin-orbit interaction within the initial 
(3(i) 10 (5/) 2 and intermediate (3<i) 9 (5/) 3 configurations. 

We have found that the spectral shapes as a function 
of photon energy are the same in both schemes. This 
looks unexpected, because the wave functions of the ini- 
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tial state are different between the schemes. Our numeri- 
cal finding can be proved if the crystal field energy is neg- 
ligible in the intermediate state, as given in Appendix. 

The obtained spectra show strong differences between 
U Miv and My edges, the former exhibits the broader 
and more complicated (multi-peak structure) profile than 
the latter as a function of photon energy. Also the spec- 
tral shapes are found quite different between the AFQ 
phase and the AFM phase. These differences as well as 
the photon polarization dependence may provide useful 
guidance to elucidate the symmetrical aspects of the ac- 
tual system when the RXS experiment will be performed 
in high pressure phase. 

The spectral shape calculated in the AFQ phase as a 
function of photon energy seems consistent with the ex- 
perimental one at ambient pressure. On the other hand, 
the spectra have to depend on the azimuthal angle, as 
shown in Fig. 5. Quite recently, Bcrnhoeft et al. studied 
the azimuthal angle dependence of the peak intensity on 
the (005) magnetic spot in URu 2 Si2, 50 having claimed 
that the spectra show no azimuthal angle dependence. 
Mixture of phases of the O x 2_ y 2 ordering and the O xy 
ordering may make the spectra less depending on the 
azimuthal angle. Otherwise this may rule out the AFQ 
type ordering possible HO candidate. 

Finally, URu2Si2 is considered to be located at the 
borderline between a localized and an itinerant descrip- 
tion. The CEF framework adopted in this work is to be 
seen only as a first approximation for the comprehensive 
understanding of the complex physics of this compound. 
The effort to this direction will be the next study. 
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Appendix: Properties of the Scattering Ampli- 
tude 

In this appendix, we derive general expressions of the 
matrix elements of the scattering amplitude in the E\ 
process, assuming that the intermediate state maintains 
rotational symmetry. This assumption seems reasonable 
when the CEF energy is negligible compared with the 
magnitude of the multiplct splitting in the intermediate 
state. Such conditions would be satisfied in many rare 
earth metal compounds. We show that the obtained ex- 
pressions lead us to a useful limit on the shape of the 
RXS spectra as a function of the photon energy; it de- 
pends only on the type of the order parameter. We will 
not discuss the case of the Ei transition, since any useful 
limits on the spectral shape have not been derived from 
a similar discussion applied to this case. Although a sim- 
ilar discussion had been done by Lovesey and Balcar, 51 
we treat rigorously the energy denominator appeared in 
the eq. (2.2), which is set aside in ref.51, in order to in- 
vestigate the energy dependence of the RXS spectra. 



A.l General results 

First, we rewrite the scattering tensor in eq. (2.2) as 

M aa ,(G,w) = -^^e- iG ^M aQ , ; >), (AT) 

3 

M aa >.j{u) = ^2E(u,A){i) \x a (j)\A}{A\x a >{j)\i) ), 

A 

(A-2) 

£( "- A)5 ^- |EA U) + ,r - < A ' 3 » 

Note that contributions from the sites belonging to the 
same sublattice are the same. We omit the subscript j 
in the following. The initial state can be expanded by a 
linear combination of | J, to) as 

|Vo) = 5>(to)|J,to), (A-4) 

within a ground multiplct of J with to being the eigen- 
value of J z . If the rotational invariance is preserved in 
the intermediate state, the eigenstates |A) are those of 
J and J z simultaneously. (Of course the intermediate 
state is occupied by a core hole.) Notice that for a given 
value of J, there exist more than one multiplets having 
the same J value but having the different energy. There- 
fore we introduce variable i telling from the multiplets 
with the same value of J, and write the intermediate 
state as |A) = \ J',m',i). Then eq. (A-2) is written as 

M aa ,{uj) = c*(™)c(m')M™ a m '(u), (A-5) 

m,m' 

Nj, J' 

«j'm=EE E e ^ ji ) 

J' i=l M=-J> 

x (J, m\x a \J', M, i){J', M, i\x a >\J, m'), 

(A-6) 

where £7j(w, J') represents E(oj,A) in eq. (A-2). The 
number of the multiplets having the value J is denoted 
by Nj. The selection rule for the E\ process confines the 
range of the summation over J' to J' = J, J±l. The eval- 
uation of the matrix clement of the type (J, m\x a \J' , to') 
is yielded by utilizing the Wigner-Eckart theorem for a 
vector operator. 52 We write down several results derived 
from M rn ''" (uj) in the following. 

First, M m < m '(w) ^ only when \m - m'\ < 2. It is 
obvious because of the nature of the dipole operators. 
Second, for to' = to, 

(a m (u>) ma(w) \ 
-ma(w) a m (w) , (A-7) 
b m (Lu) J 

a m {oj) = [J(J+ I) +to 2 ]Fj_i(w) 

+ [J(J +l)-m 2 ]Fj(u) 

+ [J(J + 3)+m 2 + 2}Fj +1 (uj), (A-8) 
b m {w) = (J 2 - m 2 )Fj_i(a;) + m 2 Fj{u) 

+ [J(J + 2)-m 2 + l}F J+1 {Lo), (A-9) 
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a(u>) = -(2J+l)F J - 1 (w) 

-Fj{lo) + (2J + 3)F j+1 {u), (A-10) 

Nj, 

Fr^^KJWKWJ'tf^E^J'), (All) 

i=l 

where (J||Vi||J') denotes the reduced matrix element of 
the set of irreducible tensor operator of the first rank. 
Note that the diagonal elements satisfy the relations 



a m (w) = a- m (u>), b m (ui) = b- m (u>). 



(A12) 



One remarkable feature is that the off-diagonal elements 
are just proportional to m, completely factorized by the 
same function a(oj). Another important result is that the 
trace is independent of m 

Tr[M m ' m H] = 2a m (uj) + b m (co) 

= J{2J + + J (J + l)Fj(w) 

+ [J(2J + 5) + l]F J+ iH, (A13) 

where TrX is the trace of a quantity X. 

Third, when \m' — m\ — 2, we can derive the following 
results. 



1 i 

M m ' m+2 = a'^a"(Lj) [ i -1 



1 -i 

M m+2 ' m = a'; n a"(uj) ( -i -1 





<C = v/(J-m)(J + m+ 1) 
x y/(J -m - 1)(J + to + 2), 



(A14) 



(A15) 



(A-16) 



a"(cu) = -Fj-!(w) + Fj(lu) - F j+1 (lu). (A17) 

Here, we can see the energy dependence is factorized by 
the same function a"(w) for every m. From these results, 
we can check a"_ i m+2 ) = a m holds, which leads us to the 
following relation 

M m > m±2 (uj) = M-( m± V'- m (cj), (A- 18) 

where either all the upper or all the lower signs are to be 
taken. 

Finally, for another m! ^ m case, that is, for \m'—m\ = 
1, the scattering amplitudes take forms 

a' m (uj) 

M m > m+1 (cj) = ( ia' m (u) 

V m (u) ib' m (u) 

(A19) 

b' m (lj) 

M m+1 ' m {cj) = | -ib' m (u) 

a' m (w) -ia' m (u) 

(A-20) 

Actual forms of a' m (w) and b' m (cu) can be obtained, but 
are not listed here because they are irrelevant to our 
present interest. We only notice that the matrix elements 
a' m (w) and b' m (ui) cannot separate their arguments u) and 



A. 2 Application to URu 2 Si2 

In the AFM phase, the relevant contributions to 
the RXS amplitude stem from the terms M m ' m (uj) — 
M~ m '~ rn (uj) oc ma(u>), which justifies the matrix form 
given by eq. (4.19). A detailed calculation reveals that 
the amplitude is proportional to the expectation value 

(Jzh 

In the AFQ phases, on the other hand, the relevant 
contributions to the RXS amplitude come from the terms 
M m ' m+2 {uj)±M m+2 ' m {uj) oc a'^a"(u) where a plus sign 
corresponds to the O x 2_ y 2 phase while a minus sign cor- 
responds to the O xy phase. This also verifies the scat- 
tering amplitudes are given by eqs. (4.16a) and (4.16b). 
In these ordering phases, the energy dependence of the 
scattering amplitude is factorized by a single function 
a(uj) in the AFM phase or a"(ui) in the AFQ phase. 
This is the reason why the shape of the RXS spectra de- 
pends only on the type of the order parameter. Note that 
this factorization is special for the E\ transition in the 
present model. If terms M m ' m±1 (aj) become relevant for 
some physical reason, the factorization breaks down due 
to the presence of the elements a' m (ui) and b' m (ui). We 
can see that the situation is more complicated for the E 2 
transition. 

Since the true nature of the hidden order phase is un- 
known, it might be useful to consider what happens if 
some octupolar phase is substantiated in URu2Si2- As- 
suming two components of the octupolar ordering T xyz 
and Tf, we can calculate the RXS signal from the oc- 
tupole phase. For instance, we choose the initial states 
as these ordering phases in the singlet scheme and apply 
the results obtained in this appendix. The result says no 
RXS signal is detected in the antiferrooctupolar phase. 
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